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Abstract. After briefly reviewing the results about polytropic stars in Palatini /('R)-theories, we 
first show how these results rely on the assumption of a regular function f{K). In particular, singular 
models allow to extend the parameter interval in which no singularity is formed. 

Furthermore, we show how the conformal metric can be matched smoothly in the cases where 
the original metric generates a singularity. In fact, the singularity comes from a singular conformal 
factor which is continuous though not differentiable at the stellar surface. This suggests that the 
correct metric to be considered as physical is the conformal metric. 

This is relevant because, even when matching the original metric is possible, the use of the 
conformal metric generates different stellar models. 
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1 Introduction 


Stellar models with polytropic equation of state (EoS) in Palatini /(*/?)-theories generically lead to 
a singularity which forms near the stellar surface for some parameters of the EoS [see 1]. Since no 
extreme physical situation is expected near the surface, this result has been used to rule out Palatini 
/(“/^l-theories. 

A further analysis [see 12] relates the singularity to the particular form of EoS, rather than to 
the dynamics of f{K). In particular, it is shown that the singularity forms at pressures which are 
unphysically low. It is thus argued that the singularity can be removed by a slight modification of 
the EoS. Since the EoS is always an approximation of the matter equations in a selected dynamical 
regime, a mathematical singularity which is unstable under modifications of the EoS needs a detailed 
discussion of the EoS itself, in order to be considered on a physical stance. 

Hereafter, we investigate the assumptions at the basis of the no-go theorem for polytropic stars 
in Palatini /(‘??)-theories. In particular, we show that allowing continuous non-analytical f{K) dy¬ 
namics, one can extend the parameter interval in which no singularity forms to cover the physical 
cases (for example, neutron stars), which were previously left out. The regularity assumption on the 
function f{K) hence appears to be essential for the no-go theorem. 

Palatini /(‘R)-theories are essentially bimetric theories with an original metric g and a conformal 
metric g. Previous analyses were done as if the conformal metric is an auxiliary one, while the original 
metric is endowed with the whole physical meaning. 

We have different reasons to consider alternatives. In the beginning of the ’70s Ehlers-Pirani- 
Schild (EPS) proposed an axiomatic approach to gravitational theories, in which the geometry of 
spacetime is derived from potentially observable quantities, instead of being assumed a priori [see 
3, 5, 9]). EPS framework clearly singles out a conformal structure, which is defined by light rays, 
and a free fall related projective structure. These two structures have to be compatible, in order to 
implement light cones as limit speed for free fall: this condition is called EPS-compatibility. When 
this holds the geometry of spacetime is a Weyl geometry, not necessarily a Eorentzian structure. 

Thus, it is natural to consider Palatini theories for a metric and a connection: the first related 
to light cones and distances, the second related to free fall. In addition, we need to restrict to those 
dynamics which implement the EPS-compatibility between conformal and projective structure. These 
theories are called extended theories of gravitation (ETG) and Palatini /(“R)-theories are in fact ETG. 

It can then be assumed that the connection f, and not the original metric, is related to free fall. 
By solving the equations and introducing the conformal metric g, the latter is thus endowed with the 
physical meaning of determining free fall of test particles. We show that in stellar models matching 
conformal metric is as easy as in standard GR and singularities at the surface do not generically form. 
It is in view of this result that one can consider the conformal metric as a candidate for the physical 
metric, even ignoring EPS framework. 

In conclusion, the matching of the conformal metric clearly spots the origin of surface singu¬ 
larities in the derivative of the conformal factor. As former analyses of no-go theorem demonstrated, 
these singularities are generic with respect to the function fiK). Moreover, matching the original or 
the conformal metric is always different, even for the parameter values which allow regular matching 
in both cases. 

1.1 No-go theorem 

In [1] it was shown that, for generic choices of analytical fCR), there are commonly used equations 
of state for which no satisfactory physical solution of the field equations can be found (apart from the 
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special case of GR), casting doubt on whether Palatini f{'R) gravity can be considered as a physically 
sound alternative to GR. 

We consider the common action for Palatini theories on a spacetime At of dimension m = 
dim(Al) given as 

s [g, r, «A] - ^ J d'”x V=g l/m) + £M(g, <A)] , (1-1) 

where we use relativistic units, i.e. c = fi - 1. Here, = SttG, G being the gravitational constant, 
g is the determinant of the metric of signature (m - 1,1), “R [g, F] = g^'^R^iv is the Ricci scalar 
built with an independent torsion-less connection f, ij/ denotes a collection of matter fields, -Cm is 
the matter Lagrangian (here assumed to be independent of f) and f{K) is an analytic function of the 
scalar curvature. 

The matter Lagrangian is assumed not to depend on the connection f . This is a standard as¬ 
sumption in Palatini f{K) gravity and it makes easier to discuss field equations hereafter. If this 
assumption is relaxed one has an extra stress tensor on the right hand side of the second field equa- 
fion (1.2b) below. These models have been discussed in fhe liferafure; [see 11, 13, 14]. However, 
generically, fhis exfra coupling prevenfs fhe connecfion f from being mefric and fhe spacetime fo 
be endowed wifh an infegrable Weyl geomefry; [see 8]. These more general models can be sfudied 
even if fhey are affecfed by holonomy problems as fhe original Weyl proposal for unified fheories of 
gravify and elecfromagnefism. 

Independenf variafions wifh respecf fo fhe mefric g^y and fhe connecfion f give fhe following 
field equafions 

f{K) R(^y) - i/(R)gpy - , (1.2a) 

. 1 

V,[V=g/'(V’'] -0, (1.2b) 


where /'(R) = d/(R)/dR, is fhe covarianf derivafive wifh respecf fo fhe independenf connecfion T 
and Tpy is fhe usual maffer sfress-energy tensor. We assume fhe latter in fhe sfandard form of a perfecf 
fluid sfress tensor, 

Tl,V^ip + P)Uf,Uy+ pgf,y, (1.3) 

where u is fhe m-velocify, p fhe energy densify and p fhe pressure of fhe fluid. The lengfh of u has fo 
be fixed so fhaf \u^ = -1. 

Equation (1.2a) can be fraced by g^'' fo obfain fhe so-called master equation 


nR)R-'^m)^K^T, ( 1 . 4 ) 

where we sef T = g^^T^y. For a generic algebraic function / fhe master equation is as well an 
algebraic equation for R, wifh T considered as a paramefer, and one can solve equation (1.4) for R fo 
obfain R = R(r). Wifh fhe latter, we can redefine / = /[R(r)] and analogously for ifs derivatives. 
The derivafive of R wifh respecf fo T is given as 


dR(r) _ __ 

dr “+ _!)//■ 


(1.5) 


Equation (1.2b) can be solved by defining a new conformal mefric g^y = ^g^y, wifh a conformal 

2 

factor (p = Consequenfly if becomes 



(1.6) 
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and can be solved with respect to the connection to obtain 



(1.7) 


In the end, by substituting equation (1.4) and the connection (1.7) into the first field equation 
(1.2a), one gets 

^ , ( 1 . 8 ) 

where G^y=R^y - ^Rg/jy is the Einstein tensor for the conformal metric g. Hereafter, tildas refer to 
quantities depending on the conformal metric g. 

In view of conformal equivalence, one can decide to write field equations in terms of g or g. 
The field equation (1.8) written in terms of g reads as 

Gf,y - J<\y , (1.9) 


where we set 


K Tf,y = 


- 2 g/iv) 


f 


m — I 1 


f 


2 (/')2 


V^/'V,/' - -VpfVPfg,, 


( 1 . 10 ) 


Here, is the covariant derivative with respect to g and □ = is the box operator of g. 

Setting the dimension m = 4, one recovers the field equations in [1]. 

It is well known that Palatini /('R)-models in vacuum are equivalent to standard Einstein models 
with a cosmological constant dictated by the master equation [see 2]. Hence, vacuum solutions are 
defined and we can assume exterior solution to be determined. To calculate the inner solution, we 
have to solve (1.9) by assuming a static spherically symmetric ansatz 


g = -A{r) dt^ + B{r) dr^ + dO^ , 


( 1 . 11 ) 


where dQ^ is the line element along the 5''”“^ angular directions. The m-velocity is then u = 'sfAjA Bq. 

We present the results of [1] for dimension m = 4 in the Appendix. We checked that the result 
is similar in dimension m = 3. 

Eurthermore, if one assumes a polytropic Equation of State (pEoS) and tries to match the inner 
and exterior (vacuum) solutions, there are values of the polytropic parameters for which singularities 
develop on the surface of the star, namely there is no matching. This result holds for generic 
analytical functions /{R.) and for physically sound values of the polytropic parameters. 


1.2 Re-examination of poly tropic star models 

In [12], the situation described in Section 1.1 was reviewed in the particular case of 

f{R) = R+^'R\ ( 1 . 12 ) 


There is a natural scale for the constant e which is of order the squared Planck length, for dimensional 
reasons. 

Under these assumptions, we can compute the divergent contributions to the scalar curvature R 
of the metric g, as p tends to 0, namely 


R ~ 


(2 - y)x2e 1 . 

(yK)^ 


+ ^^0- 


(1.13) 
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We refer to the Appendix for details and notation. Here in particular, Hq denotes the outer scalar 
curvature, Mo and vq the mass and radius of the star, K and y are the parameters of the polytropic 
equation of state (pEoS); see equation (5.8) in the Appendix. 

The scalar curvature R, in fact, diverges as p tends to 0, for 3/2 < y < 2. In common astrophys- 
ical models, for example neutron stars for which typically y = 5/3, the divergence is driven by the 
a-dimensional factor and defines a critical density scale at which the singularity is forming, namely 

1 

3-2r 

(1.14) 



For realistic neutron stars one obtains pc - 10“^'*^ [g/cm^], which is well below any reasonable 
physical density. Therefore, the existence of a curvature singularity at such extremely low densities 
should be considered as unphysical and an artifact of the idealized situation. To have an idea of 
this density scale, one single electron in the volume of a sphere with the radius of the Neptune orbit 
amounts to a density of p^ - 10“^^ [g/cm^]. 

In conclusion, the results of [12] have confirmed that the surface singularity exists, but its phys¬ 
ical feature is conjectured to depend on the details of the model. In Planck scale models, the singular¬ 
ity develops in conditions well beyond the validity of assumptions and approximations of the model 
itself. 


1.3 Problems and contents 

Generally, Palatini fCR) models together with pEoS predict a singularity at very low pressure (or den¬ 
sity), for some intervals of polytropic parameter y. Although this prevents us from modelling some 
types of stars as usually done in standard GR, one has to consider that pEoS is not a fundamental prin¬ 
ciple of physics, being instead an approximation. For example, one could assume some mechanism 
which modifies the EoS near the surface preventing the formation of the singularity. Alternatively, 
one could relax the assumption of a vacuum solution in the exterior, as we do in the following. 

Another interesting issue is related to the choice of the conformal frame. Both the primary 
no-go theorem analysis [see 1] and its re-examination [see 12] are concerned with the original metric 
g, entering the variational principle. We know, however that another natural metric g arises while 
solving the field equations. In the case of a non-regular / (or multiple zeros), the new metric g can 
gain/lose regularity, by means of the conformal factor. 

This paper is organized as follows. In Section 2, we discuss the results revised in the Appendix. 
We show that if one relaxes the assumptions about the analyticity of the function /, then the quantity 
driving the singularity development may improve -quite drastically- its behaviour. Ideally, after 
showing that the proof of the no-go theorem is not valid in these models, one should propose a 
model in which singularities can be avoided in the parameter range y € (0,2). Future investigations 
will be devoted to this topic. Section 3 examines the case of a polytropic star by embedding the 
spherical inner object into a Friedmann-Eemaitre-Robertson-Walker (FERW) background: the so- 
called Swiss-cheese model. However, the fact that singularities form at an unphysically low density 
is not enough to casually relate low pressures and singularities. Since it is difficult to obtain general 
results on the matching, one should be able to give an example, avoiding low pressures, in which 
the singularity does not appear. This turns out to be non-trivial, though investigating in this direction 
clarifies that the matching on a hypersurface has to be performed using Darmois junction conditions, 
instead of matching the metric coefficients only [see 4]. Darmois conditions are powerful because we 
can avoid assuming a coordinate system around the matching surface. In addition to this, one does 
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not have to assume that the chart is precisely the one used to express the ansatze. More specifically, 
Darmois junction conditions require that the first fundamental forms (intrinsic metrics) and the second 
fundamental forms (extrinsic curvatures), calculated in terms of the respective coordinates (outside 
and inside the matching surface), must be identical on both sides. These junction conditions allow 
us to use different coordinate systems on both sides of the matching hypersurface. In Section 4, we 
discuss how the matching procedure behaves under conformal transformations. This is motivated by 
the fact that the model is known to be similar to standard GR, with respect to g, and therefore being 
the matching possible. We find thaf fhe singularities are driven by a discontinuity in the derivative of 
the conformal factor. In the end, we draw our conclusions and physical interpretation in Section 5. 


2 Non-analytic models 


We refer hereafter to the Appendix for the detailed calculations and notation. There we review the 
original argument and find that one can define a quantity 

df 

C = (fi + p)^ ( 2 . 1 ) 

dp 

Here /' denotes the derivative of the function f{K) with respect to its argument K then evaluated at 
'R(r) with T = 3p - p. Accordingly, it is a function of (p,p). By using EoS one has p = p{p) and 
C becomes a function of p only, which is in turn a function p{r) of the coordinate r when evaluated 
along a solution. Let us then denote by Cr the derivative 


dC dp 
dp dr 


( 2 . 2 ) 


We focus now on the relevant quantity which drives the divergence of the scalar curvature of g, 
see equation (5.34) in the Appendix, namely C,.//'. Note that if one allows / to be non-analytic, then 
/' might diverge, implying consequently the convergence of the scalar curvature. 

For example, if we consider models with f{K) = c\W, then the derivative with respect to H is 
given by f = and the master equation (for m = A) reads as 


fn-2f = ci{v-2)n'' = T . 


(2.3) 


From the latter, we obtain: 


'R = 


nl/y 


, f' = 


J/y 
1 

1 y . 


(y-2)^ 


ci(y-2) 

By assuming a pEoS, we can infer the following behaviours of f with respect to p\ 


f 


ip '’' for 0 < 7 < 1 , 
\pl^ for 1 < 7 < 2 . 


From equation (2.5), it follows that 


r 


p'y for 0 < 7 < 1 , 

i-2v(y-l) , 

' p yy for 1 < 7 < I , 

l-2Ky-l) , 

p yy for j < 7 < 2 . 


(2.4) 


(2.5) 


(2.6) 
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These are non-divergent when 0 < y < 1 as well as when 1 < y < 2 and 0 < v < l/(2y - 2). 
Accordingly, if 0 < y < 1/2, these quantities are regular for a generic y € (0,2). Note that in these 
models, / is not analytic at “R = 0. This example is just to show that if one relaxes the hypothesis 
of analyticity, it is possible to regularly match the inside and outside solutions, for a value of the 
polytropic parameter y € (0,2). 


3 Swiss-Cheese: matching with FLRW 


When unphysical situations are pinpointed [see 12], one could try to build a mathematical model 
in which they do not appear. We examine such a mathematical representation, where the density p 
could not reach zero, as its minimum value is given by the mean density of the Universe, namely 
Pu » 10“^° [g/cm^]. The simplest cosmological model of spherically symmetric gravitationally 
bound systems, embedded in an expanding Universe, is the GR regime described by Einstein-Straus 
model of Universe, also called the Swiss-cheese model [see 6] . Here, a hole is created by removing 
the background material inside a spherical boundary and assuming all the mass is concentrated in 
the centre of the sphere. Mathematically, this representation is composed of a Schwarzschild metric 
smoothly connected to a Friedmann Universe. This translates into the matching across a spherical 
boundary of a FLRW metric (as exterior solution) to a Schwarzschild metric (as interior solution). 
The spherical boundary stays at a fixed coordinate radius in the FLRW frame, but it changes with 
time in the Schwarzschild system [see 4]. 

We now consider two different coordinate systems outside and inside: [r, cr, 6, 0] and \t, r, 0, (p], re¬ 
spectively. 

The general FLRW metric for the external solution can be written in spherical coordinates as 


9 9 

gF - -dr + a (r) 


dcr^ 


+ cr^{d0^ + sin^ 6d(f>^) , 


(3.1) 


where a(T) is the scale factor and K = 0, ± 1 is the constant curvature of space. 

By imposing the Darmois junction conditions on a spherical hypersurface E, we try to match smoothly 
the metric (3.1) to our ansatz for the inner region, which replaces the inner Schwarzschild metric in 
the standard Swiss-cheese model. 

The matching hypersurface E in the outer region will have equation cr - (Tq and will be parametrized 
by coordinates - (u, 6, cp). The canonical embedding of E in the outer region is given by 


i : 1, M : (u, 6, (p) (Tq, 6, cp) . 


(3.2) 


The same hypersurface E is also embedded in the inner region as 

j : H M : (u, 9, (p) {t{u), r{u), 9, cp) . (3.3) 

By standard argument, one can define fhe firsf and second fundamenfal forms on fhe hypersurface E 
wifh respecf fo each of fhe embeddings. The firsf and second fundamenfal forms wifh respecf fo fhe 
embedding i in fhe oufer region will be denoted by 

I'" - [I‘’]ABdk^ ® dk^ , ir - [ir]ABdk^ ® dk^ , (3.4) 

respectively, while fhe firsf and second fundamenfal forms wifh respecf fo fhe embedding j in fhe 
inner region will be denoted by 

r = [f]ABdk^ ® dk^ , ir' = [If ]ABdk^ ® dk^ , (3.5) 
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respectively. 

The non-null components of the two fundamental forms of S given in an implicit form, i.e. 
cr = ctq, with respect to the external space-time, are 


'[Fill = -(T'f , 
[F]22 = 0 ^ 0 - 1 , 

[F]33 = ct-q sin^e , 


[IFlii -0, 

[IF]22 ^ a CTO ^ , 

[IF] 33 ^ a cro^ sin^0 , 


(3.6) 


where we have set ^ = y 1 - Ka^ + 0. Let us remark that both the first and second fundamental 
forms are diagonal. 

As regards the inner solution, instead, we calculate the fundamental forms from the line element 
(1.11). In this case, the hypersurface S is given by f = t{u), r - r{u), and we set A = A \r{uj\, 
B - B \r{uy\. For the non-zero components we obtain 


[Fill 

mil - , 

[F]33 =r2sin2 0. 


(3.7) 


[IF'lii = -\{t' [<1)2 - {r'fB\)A' + t'ir'fAB' + 2AB{t'r” - , 


[IF ]22 = rt' yf (F-i , 

[IF] 33 = rt' sin^ 6 


(3.8) 


where we have set <1)^ = {t')^A - {r')^B. Here, A' and B' are the derivatives with respect to the 
coordinate r (then evaluated at r = r{u)), while r' and t' {r" and t") denote derivatives with respect to 
the parameter u, along the cutting hypersurface. 

Thus, we examine the matching conditions between the FLRW metric of equation (3.1) outside and 
the metric of equation (1.11) inside. The condition [F],^ = [F],y implies 


\t'fA-{r'fB = {T'f, 

a^CTg - 


(3.9) 


On the other hand, the condition [IF ]22 = [IF ]22 translates into the following equation 


ao-Q^ ^ rt' ' 


(3.10) 


This, together with equations (3.9), can be recast as 


{t'f = 


A(1 - ^B) 


{rf , 


{T'f = 


B 


(1 - 


(/)^ . 


(3.11) 


Accordingly, one cannot choose r - const since it would imply r' = 0 and t' - t' = 0, i.e. the 
matching surface Z would degenerate to a 2-surface. By using the expression for {t'f, we obtain 


{t'fr" - r't't" = ^^B 


A{2 - f^B)B' -B{1- f^B)A' 
2 A 2(1 - ^ 2^)2 


(r'4) . 


(3.12) 










Finally, we need to check [IF]ii = [II'ln, which can be written as 


2 

^ ^ A(1-B^2)2 


- 0 . 


(3.13) 


Since one cannot take / = 0, the other possibility is to have {AB)' = 0. Under this assumption, 
equation (3.13) is identically satisfied and the matching conditions read as 


\t'? = 
< (t'f = 


B 


(1-Ub) 

a(i-Ub) 

J2. 


(F)2 

(r'f 


(3.14) 


The original Swiss cheese model is recovered by setting the inner solution as Schwarzschild (and conse¬ 
quently AB - 1) and fixing t' = 1 as ansatz. Under these assumptions, one obtains 


(f')2= f"'-' 


2M 


2 2 


(r-2M)2 

,.2 


(3.15) 


The first two equations reveal us information about the cutting procedure for any provided initial conditions 
(fo, ro). The third condition, instead, constrains (fo, ro) and has to be preserved during the evolution. In fact, by 
differentiation one has 


acr 0 


Consequently, having 


one gets 


, , 2M 

^2 - 1 -H - 


= - and 


+ K 


2M 


- P ■ 


M 


2a 
— H- 
a 


a^ + K 


= 0 = p 


(3.16) 


(3.17) 


Along the matching surface the pressure is continuous and zero, while the density is positive. Hence, the 
Schwarzschild metric can be matched with a FLRW solution. 

It is quite straightforward to recover the usual polytropic matching in standard GR. In this case 
one matches an internal (non-vacuum) solution of Einstein equations with an external Schwarzschild 
solution. The matching surface is chosen to be r = ro and the matching of the second component 
of the first fundamental form implies that the radial coordinate is the same in both the interior and 
the exterior. The time coordinate t is also set to be the same inside and outside. Thus, the matching 
conditions read as: 

(a('4)'= (a«)' , 


which are verified for standard GR (see Appendix), in the special case of /' = 1, f" - 0. No 
divergence arises, since in GR C - \ and, since f - \, Cr = 0. \n general, there are two equations 
for determining the surfaces to cut and match. Initial conditions are given so that a^cr^ = at 
u = 0. Therefore, the matching is continuous if and only if the previous condition is preserved 
along evolution of [t{u), r{u)] and (AB)' = 0 along the matching surface. The real difficulty is the 
matching of the pressure on the surface, namely the matching of the p{u) value given by the exterior 
model with the inner value.The exterior dynamics depends on the function fCR) and on an external 
EoS. The function A(r) and B{r) are determined by the internal equations and one can require initial 
conditions such that AB = const. Whether this is preserved by the inner dynamics and is valid for all 
u is not clear in general. Eurther investigations should be devoted to study if it is possible to make a 
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toy model in which the matching can be done and singularities are avoided. However, it is evident 
that the matching should be done in general by connecting the first and second fundamental forms. 
In [1] it has been shown that there is no smooth matching if one assumes the surface E to be in the 
form r = vq, which is however a subclass of possible matchings. 


4 Conformal transformations 


We consider now the behaviour of a polytropic matching under conformal transformations. We decide 
to use g which obeys Einstein-like field equations, differently from g. In [10] it was shown that for a 
fluid energy momentum tensor there is a fluid effective energy-momentum tensor 

= {P + P)U^Uy+ Pg^y , (4.1) 


where the new effective pressure and density are given by 

„ _ p n{T)(P + k^T ~ _ P _ f<(T)(t> + k^T 
^ (p^ ^ (f>^ 

By assuming f{K) = 1i+ (as in Section 1.2), we obtain 

Ap + €k^T^ ^ _ Ap- eK^T^ 

^ ” 4(1 - ek^T)^ ’ ^ ” 4(1 - e/c2r)2 ■ 


(4.2) 


(4.3) 


If the original metric g is spherically symmetric, then one expects the conformal factor 0 to de¬ 
pend only on the radial coordinate, the conformal metric g being also spherically symmetric. Namely, 
we have 

g = -A{r)dt^ + B{r)dr^ + r^dO?' g - -A(f)dt^ + B{f)df^ + r^dQ^ , (4.4) 

where we have rescaled the radial coordinate = (pr^ and set 


A = (pA, 


B = 


A(p^f 


3^4 


{(p'r^ + 2(pr)^ 


B . 


(4.5) 


If we now introduce a pEoS on (p, p), we would recover the GR case, where the matching of an inner 
solution with a Schwarzschild external solution g is possible. On the other hand, if we impose a pEoS 
on (p, p) (as done in proving the no-go theorem), we can simplify for small p as follows 


|('y-l)p for 0 < y < 1 , 
yKpy for 1 < y < 2 . 


In this case, the conformal factor reads as 


cp^f'{K{T))-l-eK\2p-p)- 


j 1 - eK^(3y - 4)p 
|l - eK^p 


for 0 < y < 1 , 
for 1 < y < 2 . 


(4.6) 


(4.7) 


Eor 0 < y < 1, hence T ~ (3y - 4)p and, in view of equation (4.6), for small p the following hold 


4p -t e/f2(3y - 4)2p2 ^ 4(y _ i)p _ _ 4)2^2 

4(1 - €K^{3y - 4)p)2 ^ ^ 4(1 - €K^(3y - A)p)^ 


(4.8) 
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implying that p ~ {y - \)p. 

For 1 < y < 2, instead, T ~ 3Kp^ - p - p and for small p we can simplify as follows 


4p + eK^p^ ^ 4Kp~>' - €K^p^ 

^ " 4(1 + 6/r2p)2 ~ ^ ^ " 4(1 + eK^p)^ ~ 


(4.9) 


implying that p ~ Kp^. 

Note that an exact pEoS for (p, p) implies a varying pEoS for (p, p). However, pEoS is preserved by 
conformal transformations, at small p. To summarize, conformal transformations map Palatini /(*/?)- 
theories for g into an effective standard GR for g. The spherically symmetric ansatz, the polytropic 
equation of state and the form of source energy-momentum tensor are preserved. In the standard 
formulation for g, polytropic stars are easily produced, thus it should be possible to map g back to 
the original framework for g. However, the no-go theorem shows that there might be no matching 
available for g. The only possible explanation for this is that the conformal transformation itself is 
singular at the matching surface. The conformal factor in the interior reads as 


|l - €K^{3y - 4)p for 0 < y < 1 , 
11 - €K^p for 1 < y < 2 , 


(4.10) 


while 0 = 1 in the exterior. The conformal factor is thence continuous, though not differentiable at 
the surface. This means that a matching for g is mapped into a matching for g. In fact, the 
second fundamental forms for the matching surface with respect to g and g are related by 

[if \ In 0 7^7^ , (4.11) 

where 7^ denotes the Jacobian of the embedding of the matching surface into M. One can see that 
the transformation law depends on tp' which is discontinues at the matching surface. Eet us denote by 
[if] and [IF] the matrices with entries [if ], 4 b and [IF]^g, respectively, then the difference between 
the two matches is measured as 


if - 


1 t'(p' 


0 

, 0 


0 

-fA 

0 


0 

0 

-f A sin^(0). 


(4.12) 


which cannot vanish along the matching surface. To conclude, even when both the matchings are 
possible, matching C' with respect to g is not the same as matching C' with respect to g. The choice 
of the correct metric to be matched is a matter of interpretation of the model and its physical meaning. 


5 Conclusion and perspectives 

We showed that in Palatini /(“R)-theories problems might arise in matching inner and exterior so¬ 
lutions for stars, when using a pEoS and for analytic /(“R). In particular, for some values of the 
polytropic index (3/2 < y < 2) and a generic analytic fCR), a singularity develops near the surface 
of the star. We also proved that matching the conformal metric g is possible as in standard GR. In 
this case, the singularity is in the conformal transformation, which is instead of C^, affecting the 
matching of the second fundamental forms. 

Naturally, the interpretation of the theory gives a hint on which is the physical metric involved. If one 
assumes the free fall of particles being governed by g, then the theory is equivalent to Brans-Dicke 
theory with a potential and m = -3/2 [see 7]. However, this value is ruled out by Solar System 
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tests (though assuming no potential). In order to proceed under these assumptions, the effect of the 
potential has to be also considered. 

On the other hand, we believe that there exist good motivations for assuming that the free fall of par¬ 
ticles is ruled by g. In fact, by universality theorem [see 2] vacuum solutions are described by metrics 
which are also solutions of GR with a (suitable) cosmological constant depending on the function 
fCR). If we are modelling for example the Solar System, then we have one of these solutions and, 
since we know that the cosmological constant has no effect at these scales, we can reasonably assume 
that fiK) should be such that it corresponds to a cosmological constant which is small enough to agree 
with Solar System experiments. There are plenty of models for which the cosmological constant has 
no effect at small scales though it becomes relevant at bigger scales. 

Beside this motivation, Ehlers-Pirani-Schild (EPS) proposed an interpretation of gravitational 
theories based on potentially observable quantities (worldlines of massive particles and light rays). 
They showed that the geometry of space-time is described by a conformal class of metrics [g] and a 
projective class of connections [f]. In this framework, the free fall is associated by construction to 
the connection T and it suggests a Palatini approach [see 3, 5, 10]. Moreover, Palatini /('R)-theories 
are compatible to EPS formalism [see 9]. Within the EPS interpretation of Palatini /('R)-theories, it is 
natural to search for matching of g rather than matching of g: g is in fact the physical metric because 
it describes the free fall. In particular, this is worth to use when polytropic stars are considered. In 
fact, even when a matching of g is possible, the two matchings are not equivalent and we suggest the 
matching of g is physically more meaningful. 

Eet us remark that we did not imposed anywhere energy conditions, neither on T^y or on T^y, 
and we did not declared what is the physical meaning of the original metric g^y. The reason not to do 
that is that nothing in this paper depends on these assumptions. However, in view of EPS framework 
one should assume that the metric g is related to operational definition of distances (since choosing 
a representative in the conformal structure is in 1-to-l correspondence with choosing a definition 
for distances and time lapses). One should also investigate if ordinary quantum mechanics is the one 
for inertial observers with respect to g or g. The differences between the two metrics are probably 
too tiny to be tested experimentally and consequently also this assumption is an independent choice, 
which does not affect the content of this paper. It is our opinion that, since our protocols for distances 
strongly rely on quantum mechanics, the frame for quantum mechanics should be selected to be g as 
well. 

However, what is relevant here is that that opens a new perspective to be explored in the long 
standing issue of which is the physical frame. In the past studies, it was often, if not always, assumed 
that in the end one of the frames would have eventually emerged as the physical frame and one would 
have done everything with one metric (free fall, distances, quantum mechanics, minimal coupling, 
energy conditions, ...) while the other metric(s) would appear as auxiliary object(s) with no specific 
physical meaning. Eef us sfress fhaf fhere is anofher opfion, and if needs fo be invesfigafed, namely fhe 
possibilify fhaf differenf frames appear in fhe fheory since fhey have differenf meaning and one could 
have a reasonable model in which free fall, energy condifions, minimal coupling, quanfum mechanics 
are done in differenf frames. This possibilify corresponds fo effecfs which would be characferisfic of 
extended gravify fheories which have been almost entirely overviewed. This would give different 
frames a new and different physical meaning. 

In future investigation, one should also better examine whether it is possible to build a model 
by matching the inner solution to a EERW solution, where the density does not approach zero and no 
singularities appear, as suggested by [12]. Eurthermore, one could show whether non-analytic models 
can avoid singularities, as suggested by our analysis. 


- 12 - 


Appendix A: Inner solution in dimension m = 4 


We hereafter examine Einstein field equations for the inner solution, in dimension m - A [see 1]. The 
independent equations are obtained by plugging the metric ansatz in equation (1.11) into equation 
(1.9) and read as follows: 


B’r + B^-B 1 / f\ A 


r^B^ 


drrf (2 B'\drf 3 / d^f 
f IB) f 4\ f 


■/\2 


A'r-AB + A 1 / /, 

- ^ —Bp- -Bin- — - 
f ^ 2 \ /') 


r^A f " 2 ^ 

lA'AB - 2B'A^ + 2rA”AB - rB{A'f - rAA'B' 


2A rj f 4 i f 


,,x2 


1 


- —r^p--r^\n-—\ - 

f> ^ i \ f' I g 


f 


4A2b2 

/ 


r = 


f 


(1 B' A'\drf 3ldrf 

f 


r 2B ' 2Aj f 4\ f 


,\2 


(5.1) 


Here, primes on A,B and / denote derivatives with respect to their argument, namely r and n, respec¬ 
tively. Note that there is another equation associated to the angular coordinate 0, which is equivalent 
to the equation associated to the coordinate 0. The first two field equations can be recast in the 
following form: 


A' 

'a 

B' 

'b 


1 l\-B B 
I + (T 
1 ll-B 


2 

r f r I 


I + cr 


B 


a + P 


+ —K rp + 
r f r 


(5.2) 


where we set 


a = 




2drf 

rf 



13 = ? 


drrf 

f 





CT = 


rdrf 

2f' 


(5.3) 


The conservation law for the energy-momentum tensor, i.e. VyT^'' = 0, which is a consequence of 
Bianchi identities in equation (1.9), can be expressed as 


-f ip + P) 

1 -I- CT r[r - 2M(r)] 


M{r) + 


K^r^p 


a 

2 


[r - 2M{r)] 


(5.4) 


where we set pr=p'{r) for the derivative of the pressure with respect to the coordinate r and 


M(r)= 


r{B-\) 

2B 


(5.5) 


from which it follows that 


B = 


r 

r — 2M{r) 


(5.6) 


By fixing an EoS p - p{p) for the matter fluid, one has three unknown functions of r, namely 
{p,A,B). Equations (5.4) and (5.2) are an ODS in normal form, which in principle can be solved. 
Although the system is generally determined, the solution is in practice hard to find in an explicit 
form. In view of (5.2), we transform equation (5.5) into a differential equation for M(r), as 


M'(r) = 


1 

1 -I- CT 


K^r^p a + B 

- - H- 

2f' 2 


M{r) 

r 


{a+ (3- cr) 


(5.7) 
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It follows that the third field equation in (5.1) is identically satisfied. 


If we now consider a pEoS for the fluid 




y- 1 ’ 


fhen it follows that 


yK~y 


dp _ 

+ ^K-r 


(5.8) 


(5.9) 


Here y and K are two real constants which depend on the fluid: in particular y e (0,1) U (1,2) and 
.S' > 0. Notice that T = 3p - p and 


JO J _ 


.2 rff 




dp rn-f 


3yK~y 


— . yKy 


- 1 


(5.10) 


This stays finite while approaching the stellar surface (as p tends to zero) and has a limit which 
depends on the parameter y, namely 


f' 

Jp 


(3y - 4) for 0 < y < 1 , 

^2 r// 

1 < r < 2. 


(5.11) 


Note that for a generic analytic /, this quantity is finite. In fact, the denominator is computed at a 
zero of the vacuum master equation and it is non-vanishing if the zero is simple, which is generally 
true. Note also that in standard GR f" = 0 and = 0, for any value of polytropic parameters. 

f' 

Furthermore, we can compute the quantity fpp=-^ as 


f = 

Jpp 


.4 pm 




{f'n-f'f \ dp 


3^-11 - 


(dp 


\ dp 


3^-1 - 


,2 ftt 


3K^r 


(i-r) 


( Izz 

y\p ^ ^ 


dp\ 

d^] ^ ' 


(5.12) 


The limit of to the stellar surface as p tends to zero depends again on the value of y 


f 

Jp, 


pp 


(- (3y-4)2 

-iK-f" d-yf 


fn-f 

3^V" 




for 0 < y < I , 
for 5 < y < 1 , 


2 r-2 


y{l-y)Kyp y forl<y<2 
though it does not stay always finite in this case. It follows that: 


(5.13) 


ip + p)-r = 

dp 


( i 

py 

T 

K~y 


-I- 


yp 
y- 1 


1 py 


-1 


-I- 


y K-y y-^ 


2_1 
1 py ^ 


2py 


+ 


yp 


y K~y {y-l)Ky (T - 1)^ 


(5.14) 


which approaches zero as p tends to zero, if 0 < y < 2. 
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We multiply equation (5.4) by fp=df'/dp and use equation (5.3) to get a quadratic expression 
in drf ■ Solutions of the latter are 


„ _ -drf'iC - f')(r - 2M) ± 
~ r2(3C - 4/')(r - 2M) 


(5.15) 


where we set 


df ! dp\ dp 


K^fr\r - 2M) {sf(C - f)\r - 2M) - C{4f - 3C) pK^p - f"R + f)? + 4/'m]) . 


(5.16) 


To summarize, when the stellar surface is approached from the inside (i.e. as r tends to ), then p 
tends to zero, so that also p tends to zero. This implies that T approaches zero as well, while / tends 
to /'‘Ro/ 2, where Hq is a solution of master equation in vacuum. Hence, C approaches zero, with 
df'/dp being finite. 

We set Mo=M{ro) so that, in view of equation (5.5), we can write B\rg = (1 - 2Mo/ro)~^. It follows 
that 

2A|,„=16r2(ro - 2Mo)\f'f . (5.17) 

The quantity C is a function of p (since the scalar curvature 7? is a function of T = 3p - p because of 
the master equation and the pressure density p can be eliminated by using the EoS). Moreover, since 
the pressure is a function of r on-shell, one can define the derivative of C with respect to r, which will 
be denoted by Cr- Thus, by setting Cr = prdC/dp, equation (5.15) gives either 


drflro ^ , drrf'lro = 


ro 


i4-'Rorl)Cr ^ 2f' 
8 (ro - 2Mo) 


+ 


(5.18) 


or 


drflro = 0 , 

In the case of solution (5.19), one also has 


_{'Rorl-SMo)Cr 

^rrj Iro — 


8ro (ro - 2Mq) 


(5.19) 


A' 

~A 


'•o 


8Mo - liorl B' 

4ro(ro-2Mo)’ B 


ro 


(2f + roCr/morl - 8Mo) 
8ro/' (ro - 2Mo) 


(5.20) 


and 


= 0 , 


M' 


ro 


2f"Rorl + {Korl - 8Mo)C, 
16/' 


(5.21) 


The exterior solution g is obtained by setting p - 0 (and consequently p = 0). The master equation 
simplifies to f"R - 2/ = 0, which sets a constant value for K = Ro in a quantized set encoded by the 
function /. This in turn implies that drf = 0 and drrf - 0. Vacuum field equations are thence given 
by 

^ -Agpy , where A=^ |ro - ^ ^Ro . (5.22) 


It is then reasonable to assume for the exterior solution the well-known Schwarzschild(-AdS) solution 
with cosmological constant: 

g = -A(r)df^ -I- B(r)dr^ -I- fdQ^ , (5.23) 
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where 


1 


(5.24) 


2m Rn o 

A(r)^l- 

r 12 


and B(r)= 


A(r) 


The exterior solution evaluated at the surface, namely r = ro, gives 


A' 

2(r3Ro-12m) 

A 

ro ro (r^Ro - 12ro + 24m) ’ 

B' 

2 (r^Ro - 12m) 

B 

ro ro (r^Ro - 12ro + 24m) 


(5.25) 


The exterior and interior solution need to match on the stellar surface. This can be obtained by 
imposing continuity at r = ro of the coefficients B(ro) = B(ro) and (A'IA){ro) = (A'/A){ro), thus by 
matching the constants 


'Ro - Ro 


and 


Mo -m + 


24 


(5.26) 


As a consequence, A', p, p' and drf are also continuous. We can now consider the following 


B' 

¥ 


B' 

¥ 


rlRo 


12m Cr 


rlRo 


12ro + 24m /' 


(5.27) 


To obtain continuity, one needs either to fine tune the parameters (so that r^'Ro - 12m) or have Cr! f 
approaching zero as r tends to ro. If / is assumed to be regular at ro, that amounts to have that Cr 
tends to zero as r tends to ro. Thus, B', M' and drrf' are continuous if Cy tends to zero. Notice 
however, that more generally, if / diverges at ro, then Cy is allowed to diverge, provided that the ratio 
Cylf tends to zero. We focus now on C,. at the surface: 


dC 

Cy = —Pr = 
dp 


(p + p)fpp + I + 1 l/p 


dp 


P ^ 


f 

J p 


r i ^ 

py yp 


PP 


+ 


Ry 


+ fp 


/ IzZ 
p y 


+ 


r 


yK 


y r- 1 


where we used the pEoS. In addition, the following hold: 


Pr , (5.28) 


^ dp dp 


( IzZ 

p y ^ 1 


/¥ 


f = 

Jpp 


yRy 

Izz 

n r n y 1 

+ 2- - -r + 




y'^Ry y{y-\)Ry (T - 1)^ 


i-2y 

(1 - y)p ^ , 

Jpp \ Jp ’ 


y^Ry 


(5.29) 


where we set fpp=A^f'jdp'^ and fpp=(Mf j. By using equation (5.4), one gets: 


Cy 


Roro - 8Mo (T - 1) ciCp)// + o-2{p)fpp 


8ro(ro - 2Mo) 


yHr- 1 )"^ 


(5.30) 


where by ~ we denote equality of dominant terms, as p tends to zero. The terms cri(p) and cri(p) 
read as: 


o-i(p) ^ -(y - 2)(y - lyR yp y^ 


)H 

_4 4-2y 


9 -2 2^ 9 _1 i A 

■ y{y - 4)(7 - \) R ypy + y (y - 1)(7 + 2 )R ypy + py 


(Tiip) - (y- ^)'^R ypy'+ Ay{y - 1 )^^" ypy+ Ay^{y^ - \)R ypy+ 67^(7 - if'R ypy + p^y‘^. 

(5.31) 
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Note that, referring to equation (5.12), we obtain that fppCr 2 {p) tends to zero for generic vales of y in 
the interval (0,2). As we showed that f'p is finite, for 0 < y < 3/2 Cr tends to zero. However, for 
'ijl <y <2,Cr diverges and reads as 


^ •Rpro - 8Mo 2 -y 

8ro(ro - 2Mo) ^ ^ ' 

The Ricci scalar of the metric g at the surface is given by 

r{rA' + 4A)(BA' - AB') + lBr^{AA" - A'^) - 4A^B^ + 4BA^ 

2r2A2B2 

Furthermore, we obtain that at the surface 


R 


^ rp<Q - 8Mo Cr 


8r2 

O /-0 


/' 


+ ‘Ro 


(5.32) 


(5.33) 


(5.34) 


This means that, for 3/2 < y < 2 and regular /, Cy! f is divergent as the scalar curvature. Hence, the 
singularity at the surface is not a coordinate singularity. We checked that also in the case of dimension 
m = 3, the singularity at the surface is not a coordinate singularity. 
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